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ABSTRACT 

The determination of Li and proton-capture element abundances in globular cluster (GC) giants 
allows us to constrain several key questions on the multiple population scenarios in GCs, from forma- 
tion and early evolution, to pollution and dilution mechanisms. In this Letter we present our results 
on Li abundances for a large sample of giants in the intermediate-metallicity GC NGC 6121 (M4), for 
which Na and O have been already determined by Marino et al. The stars analyzed are both below 
and above the red giant branch bump luminosity. We found that the first and second generation 
stars share the same Li content, suggesting that a Li production must have occurred. This is a strong 
observational evidence providing support for the scenario in which asymptotic giant branch stars are 
GC polluters. 

Subject headings: globular clusters: individual (NGC 6121)- stars: abundances - stars: individual 
(Population II) 



1. INTRODUCTION 

The presence of multiple populations as a characteriz- 
ing property of globular clusters (GCs) is widely accepted 
nowadays. Since the first photometric and spectroscopic 
studies, which revealed multiple sequences (e.g. Dickens 
& Woolley, 1967; Lee et al. 1999, 2009; Bedin et al. 2004) 
and/or light element (anti) correlations (see Gratton et 
al. 2004), it became evident that GC stars are neither 
coeval nor chemically homogeneous. Hence, GCs host 
at least two stellar generations. Thanks to the advent 
of 8—10 m class telescopes, which also allowed target- 
ing fainter main-sequence (MS) stars, several abundance 
studies have shown that chemical (anti)correlations are 
also present in unevolved (MS/turnoff (TO)) or scarcely 
evolved (SGB) stars (Gratton et al. 2001; Carretta et al. 
2004; Pasquini et al. 2005). This evidence implies that a 
previous generation of stars has activated CNO, NeNa, 
and MgAl cycles in their interiors, in order to deplete 
O and Mg and enhance Na and Al, respectively. The 
origin/nature of such kind of stars is still debated with 
two coexistent scenarios: (1) intermediate-mass asymp- 
totic giant branch (AGB) stars in hot bottom burning 
phase (Ventura & D'Antona 2009) and (2) fast rotating 
massive stars (FRMS; Decressin et al. 2007). 

Our group has recently carried out an extensive survey, 
focusing on the determination of proton-capture elements 
in 19 GCs (Carretta et al. 2009c), with the main objec- 
tive of discovering and understanding the nature (and the 
extent) of the chemical (anti)correlations and their link 
with global cluster parameters (horizontal branch (HB) 
morphology, metallicity, etc.). The large sample of stars 
(~1200 GC members), analyzed in a very homogeneous 
and accurate way, revealed that while the Na— O anti- 
correlation is present in all the GCs (i.e. the second gen- 
eration is not a "perturbation"), the shape of the Na— O 
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distribution varies from cluster to cluster (Carretta et al. 
2009c). On the other hand, the analysis of UVES spec- 
tra for ~ 200 stars by Carretta et al. (2009b) has shown 
that the Mg— Al anticorrelation is not present in all GCs. 
Both these indications suggest that the typical polluter 
masses change from cluster to cluster: this variation is 
apparently driven by a combination of cluster luminosity 
and metallicity. 

In this context, lithium abundances offer a complemen- 
tary approach to p-capture elements allowing to address 
several important issues. If no Li is produced by the pol- 
luters, the multiple population scenario predicts that Li 
and O are positively correlated, while Li and Na anti- 
correlated. Na-poor, O/Li-rich stars are the first pop- 
ulation born in the cluster (they share the same chem- 
ical composition of field stars at the same metallicity), 
and Na-rich, O/Li-poor stars constitute the second gen- 
eration. Within the same hypothesis, Li is an excellent 
tracer of the dilution process acting within each star: 
only through Li abundance determinations we can deter- 
mine the amount of pristine (and of polluted) material 
present in each star. In particular, we hope to answer two 
fundamental questions: do 100% polluted stars (Li~0) 
exist or even the most extreme population still contains 
a certain fraction of primordial matter? Also, is the min- 
imum measurable Li content the same for all GCs or does 
it vary from cluster to cluster? 

On the other hand, Li offers the exciting chance to 
observationally constrain the nature of the polluters. If 
the progenitors of second generation stars are FRMS, 
they have destroyed their original Li content. On the 
other hand, if AGB stars are responsible for intra-cluster 
pollution, they may have non-negligible Li yield, given 
the Li production via the " 7 Be transport" mechanism 
(Cameron & Fowler 1971). As a consequence, we can 
reveal if AGB stars are responsible for GC pollution 
it:hrough two main observational evidences: (1) the pres- 
ence of very Li-rich stars among GC populations and 
(2) the lack of Li— Na anti-correlation (or Li— O corre- 
lation), with the second generation stars also showing a 
rather high Li content. In a recent work, we focused on 
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~90 TO stars belonging to the metal-rich GC 47 Tuc: in 
this case, likely because of the high metallicity, a large 
star-to-star scatter in Li abundances erases any Li— Na 
anticorrelation, while Li and O appear to be only weakly 
positively correlated (D'Orazi et al. 2010). The cluster 
seems to display a different behaviour from NGC 6752 
and NGC 6397: for the first one, Pasquini et al. (2005) 
detected a significant Li— Na anticorrelation (and also 
Li— O correlation and Li— N anticorrelation). Concern- 
ing the second cluster, the situation seems more complex 
since Lind et al. (2009) detected a quite constant value in 
Li abundances, with only three stars (out of 100) driving 
a hint of Li— Na anticorrelation. Enlarging the sample 
of simultaneous determinations of Li and p-capture el- 
ements in GCs is hence of paramount importance: in 
this Letter we present Li results on the intermediate- 
metallicity ([Fe/H]=-1.18; Carretta et al. 2009a) cluster 
NGC 6121 (M4), by analyzing the same sample of a hun- 
dred red giant branch (RGB) stars studied by Marino 
et al. (2008). As found by Marino and coworkers, M4 
hosts two distinct populations of stars, mainly character- 
ized by a different sodium content (i.e., the Na-rich and 
Na-poor groups) , and defining different sequences on the 
color— magnitude diagram U versus (U — B). We derived 
Li abundances for the stars belonging to the two groups, 
and located both below and above the RGB bump lumi- 
nosity; this evolutionary stage plays a fundamental role 
in this context. 

Theoretical models (Iben 1967), confirmed by obser- 
vations of field stars by Gratton et al. (2000), predict a 
depletion in Li due to the first dredge-up (1DUP) of a fac- 
tor of ~ 20 at the base of the SGB branch. On the lower 
RGB, below the bump luminosity, the molecular weight 
gradient associated with the H abundance jump acts as 
a barrier that prevents further extra mixing (Charbonnel 
et al. 1994); then the Li abundance remains constant un- 
til the RGB bump is reached. At this stage, the H shell 
reaches and cancels this discontinuity and non-standard 
mixing processes (non-convective extramixing, see, e.g., 
Charbonnel & Zahn 2007) cause the total destruction of 
the remaining Li. 

2. SAMPLE AND ANALYSIS 

Our sample consists of 104 RGB stars, whose spectra 
were acquired with FLAMES on VLT/UT2 (Pasquini et 
al. 2002) with the fiber link to the high-resolution spec- 
trograph UVES (R ~ 40,000). A detailed description 
of observations, target properties, and data reduction 
and derivation of atmospheric parameters is provided in 
Marino et al. (2008). Adopting Kurucz (1993) model at- 
mospheres and using the ROSA abundance code (Grat- 
ton, 1988), we derived Li abundances by means of a spec- 
tral synthesis of the Li I resonance doublet at 6708 A. We 
changed the CN values for the two different groups of Na- 
rich and Na-poor stars (threshold value at [Na/Fe]=0.2 
dex) in order to optimize the synthesis best fit and to ac- 
count for the CN enhancement in the Na-rich population 
(see Marino et al. 2008). 

Abundances for Nai along with stellar parameters and 
metallicity are the ones presented in Marino et al. (2008). 
Concerning Li abundance, error estimates have been 
computed in the same fashion as described in D'Orazi 
et al. (2010), taking into account both stellar parameter 
and best-fit uncertainties; for errors in Na (internal and 



systematic), we refer the reader to Marino et al. (2008). 
Stellar parameters and abundances are given in Table Q] 
(completely available in electronic version through CDS). 

3. RESULTS AND DISCUSSION 

In Figure [TJ we show the Li abundances as a function 
of the absolute magnitude My for all our sample stars: 
as expected, Li disappears above the bump luminosity 
(M v =-0.05±0.10; Ferraro et al. 1999). If we focus on 
the region below the bump level (at the left side of the 
dashed line in Figure [lj , there is no systematic differ- 
ence in Li abundances between Na-rich (filled squares) 
and Na-poor stars (empty symbols). However, when we 
look at the diagram as a whole we can see a different 
drop in the Li content with magnitude for the two popu- 
lations. Specifically, while Li seems to have a gentle de- 
crease with luminosity for the Na-poor stars, the Na-rich 
group presents a very abrupt decline, i.e., at the bump 
luminosity Li suddenly disappears. This fact, which re- 
flects different timescales for mixing and hence for Li 
depletion 3 , suggests a structural difference between Na- 
rich and Na-poor stars; however, no current theoretical 
model predicts such a behavior and we cannot provide 
a satisfactory explanation to date. In this context, we 
mention that the so-called "thermohaline mixing" has 
been proposed responsible for non-canonical mixing act- 
ing at the RGB bump (see, e.g., Eggleton et al. 2007; 
Charbonnel & Zahn 2007). Eggleton and coworkers sug- 
gested that the molecular weight inversion created by 
the 3 He( 3 He,2p) 4 He reaction could be the cause of such 
a mixing: why Na-rich and Na-poor stars should be dif- 
ferentially affected by this kind of mixing is not obvious 
and our result could be the input for further theoretical 
and observational investigations in this direction. 

By considering only the stars fainter than the bump 
luminosity, we show Li abundances as a function of Na 
in Figure [2] as one can see, there is no Li— Na anti- 
correlation, with second generation stars ([Na/Fe]>0.2 
dex) sharing the same Li content of the primordial pop- 
ulation. As an example, we show in Figure [3] the spectra, 
around the Li I region, for the two most extreme cases in 
Na abundances. The two stars, with [Na/Fe]=— 0.02 and 
[Na/Fe]=+0.43, respectively, show identical Li features 
(note that the stars have very similar parameters, and 
the same line strength reflects the same Li abundance). 
The average Li abundances are logn(Li)=1.336±0.023 
(rms 0.062) and logri(Li)=1.387±0.038 (rms 0.136), re- 
spectively for Na-poor and Na-rich stars: although we 
derive the same Li abundance for the two populations, 
it is interesting to note that Na-rich stars have a larger 
scatter in Li with respect to Na-poor ones. A one-tailed 
Fisher test returns a 5% probability such that a differ- 
ence can be obtained by chance. 

A natural explanation for this similarity in Li content 
between first and second generation stars (with the last 
showing a larger scatter) is a Li production. In fact, if 
a decrease in O of ~50%/60% occurred (as derived e.g. 

3 We estimated the e-folding time for Li abundance using the 
tracks by Bertelli et al. (2008). Na-poor stars reduce to about 
a factor of 2 their Li content in 0.17 mag , which corresponds 
to fti 10 Myr (this is indeed the time required by a 0.8 Mq star 
to become brighter by 0.17 magnitudes after it has left the RGB 
bump). For the Na-poor stars, this time is smaller by at least an 
order of magnitude, i.e. <2Myr. 
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by Marino et al. 2008 and Carretta et al. 2010), also Li 
must have been depleted. In a recent work, D'Antona 
& Ventura (2010) have presented the expected Li pro- 
duction as a function of the polluter mass (AGB stars) 
for metallicity .2=0.001. Looking at their Figure 5, one 
can see that a very low mass AGB polluter (i.e., w4 Mq) 
can produce a moderate Li content with values very close 
to the Li plateau (logro(Li)~2.2-2.3). After considering 
a depletion of a factor of ^20 at the 1DUP, this result 
agrees very well with our values (i.e. logn(Li)~l. 3-1.4). 
As also briefly explained in the Section[T](and widely dis- 
cussed in Carretta et al. 2010), there are further indica- 
tions that only low-mass polluters could have contributed 
to the observed chemical pattern in M4: (1) an almost 
"vertical" Na— O anticorrelation, with very small oxygen 
variation (depletion); and (2) the lack of Mg— Al anticor- 
relation 4 , which in fact requires high mass polluters for 
the activation of higher temperature cycles (T ~65 MK; 
Prantzos & Charbonnel 2006). A similar case could have 
occurred for NGC 6397, where according to Pasquini et 
al. (2008), two stars differ by ~0.6 dex in O, but have 
the same "normal" Li (logn(Li)=2.2). Also, in a recent 
work Lind et al. (2009), based on a sample of ^100 MS 
and early SGB/stars, found no difference in Li abun- 
dances between Na-rich and Na-poor stars with only two 
stars driving a Li— Na anticorrelation. They concluded 
that Li content is independent of intra-cluster pollution; 
however the Na— O distribution points out to a certain 
(though small) degree of oxygen depletion and, as a con- 
sequence, of Li destruction as well. Hence, if first and 
second generation stars share the same Li abundances, 
a Li production should also be required for this cluster. 
Along with a difference in metallicity of ~0.9 dex, the 
two clusters, NGC 6397 and M 4, have both quite small 
integrated magnitudes (i.e., mass) with Myt=— 6.63 and 
Mvt=— 7.20, respectively (Harris 1996). The similarity 
in masses between these two GCs also seems to suggest 
a similar "typical" polluter for both M4 and NGC 6397, 
with the requirement to have in both cases neither very 
high mass polluters (no extended MgAl/NaO anticorre- 
lations, very little He enhancement 5 , and no Li— Na anti- 
correlation) nor low mass polluters (<4 Mq), otherwise 
the C+N+O is not constant and/or s-process variations 
should be present (see Ivans 1999; Yong et al. 2008). 

The more massive GC NGC 6752 (M vt =-7.73) could 
present a different behavior. We might speculate that 
only a very low Li production from higher mass polluters 
of «5-6Af (see Figure 5 of D'Antona & Ventura, 2010) 

4 We note that Marino et al. (2008) found evidence for a small 
increase (~0.10 dex) of Al with Na, with Na-rich stars also slightly 
Mg depleted (see their Table 7). Evidence for a Na— Al correlation 
was also found by Ivans and coworkers. In any case, the Al vari- 
ation is small, and no Mg— Al anticorrelation has been observed 
among M4 stars by Ivans et al. (1999), Marino et al. (2008), and 
Carretta et al. (2009b). Also the Mg isotope ratios, as derived by 
Yong et al. (2008), show no variation within M4 with values very 
close to solar ones: the same authors concluded that this is not 
surprising due to the very little Al variation in this GC. 

5 Although He is the main product of hot H-burning through 
CNO cycle, the relationship between He and p-capture elements is 
quite complex. For "extreme" GCs, like NGC 2808 and ajCen (with 
extended NaO-MgAl anticorrelations, multiple MSs and peculiar 
HB morphology) the Y content can reach up to Ri0.40. For all 
the other GCs, even if they show the NaO anticorrelations, large 
differences in Y are not necessary (see for details Gratton et al. 
2010). 
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Fig. 1. — Li as a function of absolute magnitudes for Na-rich 
(filled squares) and Na-poor stars (empty squares). The dashed 
line marks the bump luminosity as derived by Ferraro et al. (1999), 
while the solid line is an eye fit to Na-poor population. (A color 
version of this figure is available in the online journal.) 

does not erase Li— Na anticorrelation for this cluster 6 . 
Note in fact that NGC 6752 presents an extended Na— O 
anticorrelation and a large variation in Al (i.e., the MgAl 
chain was active in the polluter stars). 

On the other hand, it seems very difficult to discrim- 
inate the nature of polluters and their properties for 47 
Tuc: maybe this GC was similar to NGC 6752 but the 
intrinsic scatter in Li abundance, independent of intra- 
cluster pollution, washes out the fossil imprint by the pre- 
vious generation of polluter stars (D'Orazi et al. 2010). 



Given the large uncertainties linked to model predictions 
(cross sections, mass loss law, overshooting, convection 
treatment, etc.), here our general aim is to provide a 
"qualitative" comparison between theoretical prescrip- 
tions by D'Antona & Ventura (2009) and observational 
evidence. Also, we stress that our result, based on only a 
few objects, needs to be confirmed by including a larger 
number of clusters and of star per cluster. However, we 
think that our data provide a quite robust observational 
evidence of AGB stars as responsible for GC pollution. 



4. SUMMARY AND CONCLUSIONS 

6 The trend of Li production with polluter mass is not linear (see 
Figure 5 of D'Antona &; Ventura, 2010): the low mass polluters 
(3.5—4.5 Mq show a Li production close to the Plateau values, 
while in the higher masses regime (5—7 Mq) the Li yields become 
smaller (a "v-like" distribution). Moving to very high-mass stars 
(7-8 Mq), the Li production reaches extremely high values, up to 
logn(Li) ks4 dex. 
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Fig. 2. — Li vs. Na for stars below the luminosity bump. Error 
bars for Na come from Marino et al. (2008); the uncertainties in 
Li are due to errors on best fit and effective temperatures. 
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Fig. 3. — Comparison of two spectra, near the Li doublet at 6708 
A, for the two extreme cases for Na abundances. Solid and dashed 
lines are for [Na/Fe] =-0.02 and [Na/Fe] =+0.43, respectively. (A 
color version of this figure is available in the online journal.) 



We report in this Letter Li abundances for a sample 
of ~100 giants, both below and above the RGB bump 
luminosity, belonging to the GC NGC 6121 (M4). The 
main purpose of our work was the study of the correla- 
tion (if any) between Li content and elements involved 
in p-capture reactions, the last ones already presented 
for the very same stars by Marino et al. (2008). The 
principal results we obtained in our investigations can 
be summarized as follows. 



1. As expected, Li tends to disappear as the stars 
reach the RGB bump luminosity; however the Na- 
rich and Na-poor stars show very different trend 
of Li with magnitude. Specifically, while the de- 
cline of Na-poor stars with My is rather smooth, 
there is a very brusque decrease for Na-rich stars. 
The so-called "thermohaline mixing" , which seems 
responsible for extramixing processes at the RGB 
bump, is not predicted to have different outcomes 
in Na-poor and Na-rich stars. Further investiga- 
tions of this aspect are mandatory both from ob- 
servational and theoretical points of view. New ob- 
servations, focusing on Li determination along the 
RGB in several GCs, are necessary to assess if M4 
is a "peculiar" case or other GCs share the same 
behavior; as a consequence model predictions could 
be revised in this sense. 



2. M4 does not show any Li— Na anticorrelation, with 
first and second generation stars having almost 
the same Li content. Along with similarities in 
Li abundances, the larger scatter found in Na-rich 
stars, indicates that a Li production, from the pre- 
vious generation of polluters, must have happened. 
This provides support to intermediate-mass AGB 
stars responsible for intra-cluster contamination, 
since FRMS can only destroy Li. 
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